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Summary 

We  have  collected  seismic,  geologic,  gravity  and  potential  field  data  to  develop  a 
model  of  the  North  African  lithosphere  via  an  integrated  approach.  We  are  currently 
modeling  the  gravity  data  along  three  transects  across  North  Africa:  1)  from  north  central 
Egypt  across  Libya  to  the  Hoggar  Uplift,  2)  from  central  Tunisia  across  northern  Algeria 
to  the  Atlantic  coast  of  Morocco,  3)  from  north  central  Algeria  across  the  Hoggar  Uplift 
to  northern  Chad.  These  transects  represent  the  major  propagation  paths  between 
earthquakes  and  nuclear  explosions  and  seismograph  stations  located  within  North 
Africa.  The  gravity  models  are  constrained  wherever  possible  by  geologic  and  other 
geophysical  information.  Crustal  thicknesses  we  have  obtained  for  the  Hoggar  Uplift 
(38  kni)  agree  well  with  preliminary  results  of  seismological  studies  by  other  research 
groups.  Our  modeling  results  also  suggest  that  the  mantle  anomaly  under  the  Hoggar 
Uplift  is  relatively  small  (200-300  km  wide,  less  than  60  km  in  vertical  extent)  compared 
to  mantle  anomalies  associated  with  other  hot  spots  and  rift  zones. 

We  have  used  our  modeling  results  in  conjunction  with  other  geophysical  studies 
to  construct  a  preliminary  map  of  crustal  thickness  for  Africa.  In  most  places  away  from 
the  continental  margin,  the  crust  in  North  Africa  is  35-40  km  thick.  Crustal  thinning  is 
associated  with  the  Sirt  Basin  and  the  West  and  Central  African  rift  system.  Moderate 
crustal  thickening  is  associated  with  the  Atlas  Mountains. 

In  addition  to  our  gravity  studies,  we  have  been  compiling  a  collection  of 
seismograms  from  Worldwide  Standardized  Seismograph  Network  (WWSSN)  stations 
located  at  regional  distances  for  magnitude  >  5.0  earthquakes  occurring  within  North 
Africa  (between  1963  and  1988)  and  their  associated  aftershocks.  We  have  selected 
these  earthquake  sequences  for  study  since  many  of  the  mainshocks  have  known  focal 
mechanisms  determined  from  teleseismic  waveform  modeling  or  first  motion  studies  of 
regional,  teleseismic  and/or  local  network  data.  Many  of  the  events  produced  numerous 
aftershocks  of  varying  magnitudes.  Using  the  lithospheric  models  obtained  from  the 
gravity  modeling  as  first  approximations  of  the  velocity  structure,  we  will  model  the 
regional  waveforms  of  these  earthquake  sequences.  We  have  also  collected  delay  time, 
travel  time,  source  parameter  and  magnitude  information  for  these  sequences  for  further 
analysis. 
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Introduction 

North  Africa  has  been  the  object  of  considerable  oil  exploration  and  a  fair  number 
of  papers  containing  geological  interpretations  of  upper  crustal  structure  have  been 
published.  However,  very  few  papers  contain  basic  data  which  are  relevant  on  the  scale 
of  a  lithospheric  investigation.  Our  task  has  been  to  develop  a  model  of  North  African 
lithosphere  via  an  integrated  analysis  of  seismic,  potential  field,  and  geologic  data.  In 
particular,  we  have  concentrated  on  constructing  detailed  2-D  models  from  known 
earthquake  source  regions  to  key  seismic  monitoring  stations  in  the  region.  An 
outgrowth  of  this  effort  is  a  data  base  of  geological  and  geophysical  information  which 
will  be  made  available  to  the  scientific  community  through  electronic  access. 

'  Although  data  are  sparse  in  many  regions,  it  is  clear  that  North  Africa  has 
experienced  a  complex  evolution  (Figure  1).  The  structure  of  the  lithosphere  presently 
must  bear  the  signature  of  these  tectonic  events,  and  seismic  waves  propagating  through 
it  must  be  affected.  The  Pan- African  orogeny  (750-500  Ma)  may  represent  the  time  at 
which  the  shield  areas  formed  by  terrane  accretion  (Black  et  al.  1994).  This  could  be 
similar  to  the  formation  of  the  China  region  in  the  Paleozoic  and  Mesozoic.  During  the 
Paleozoic,  the  Caledonian  and  Hercynian  orogenies  affected  the  region  and  formed  a 
series  of  basins  and  uplifts  in  a  manner  similar  to  the  evolution  of  the  midcontinent 
region  of  the  U.S.  during  this  same  period.  In  the  Cretaceous,  Africa  was  subjected  to 
widespread  rifting  and  nearly  split  apart  (Fairhead,  1992).  During  the  Tertiary,  the 
Alpine  orogeny  affected  North  Africa,  most  notably  forming  the  Atlas  Mountains. 
Volcanic  centers  such  as  Hoggar  and  Tibesti  formed  on  topographic  uplifts  which  are 
believed  to  related  to  mantle  plumes  (Figure  1). 

Our  approach  to  constructing  a  lithospheric  model  and  data  base  for  North  Africa 
has  been  to  integrate  information  such  as  well  logs,  gravity  measurements,  aeromagnetic 
surveys,  heat  flow  measurements,  and  geologic  maps  to  first  develop  lithospheric 
models  along  key  profiles  in  North  Africa.  We  then  will  use  these  models  as  a  starting 
point  for  modeling  seismograms  at  regional  distances,  refining  the  models  as  the  seismic 
data  require.  In  following  subsections  we  describe  our  data  collection  process.  We  then 
discuss  the  results  of  our  gravity  modeling  and  analysis.  The  appendices  describe  the 
data  bases  we  have  collected  and  how  they  can  be  accessed  electronically. 


Methods,  Assumptions  and  Procedures 

Geological  Data  Base 

We  have  digitized  the  locations  of  numerous  geological  features  within  the  North 
African  region  to  aid  in  our  interpretation  of  gravity  and  seismic  data.  These  features 
include:  surface  exposures  of  Precambrian  rocks  in  the  major  uplifts,  depth  to  Paleozoic 
and  Mesozoic  basement  within  the  Sirt  Basin  of  Libya,  surface  faults  within  the  Sirt 
Basin,  the  location  of  Pan-African  suture  zones  and  Cretaceous  and  Cenozoic  age  rifts, 
and  the  boundaries  of  West  and  Central  African  rift  basins.  Maps  of  these  features  are 
shown  in  Appendix  A,  along  with  details  on  how  these  outlines  can  be  accessed 
electronically. 
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Figure  1.  Generalized  tectonic  maps  of  northern  Africa  A-northwest;  B-northeast 


Gravity  Data  Base 

We  have  collected  gravity  data  from  a  number  of  sources.  First  we  began  with 
the  unclassified  data  from  the  global  data  base  maintained  by  the  Defense  Mapping 
Agency.  For  these  data,  we  have  complete  information  including  elevation  and  there  are 
over  36,000  data  points.  To  this  we  have  added  data  from  Libya  which  were  compiled 
by  one  of  our  graduate  students  (Suleiman,  1993).  These  data  (over  5000  points)  come 
from  a  variety  of  industrial  sources,  and  thus,  complete  information  is  lacking.  We  do 
have  latitude,  longitude,  and  Bouguer  anomaly  value  (density  -  2.67  gm/cc).  Details 
concerning  these  data  are  available  in  Appendix  B. 

A  third  and  very  important  set  of  data  was  obtained  in  the  form  of  a  5'x5' 
Bouguer  gravity  anomaly  grid  from  GETECH,  which  is  associated  with  Leeds 
University.  These  data  are  for  the  region  between  0°  and  37°N  and  12°W  and  40°E.  We 
are  only  free  to  use  these  data  in  the  form  of  maps  which  are  contoured  at  intervals  of  10 
mGal  or  greater.  Gravity  station  locations  were  also  provided  so  that  we  can  clearly  see 
the  constraints  that  these  data  provide.  In  addition,  the  station  elevations  were  used  by 
GETECH  to  correct  errors  in  the  ETOP05  topographic  data  base,  and  we  are  using  this 
commercial  product  in  our  study.  The  ETOPO  5  data  base  is  available  from  the  National 
Geophysical  Data  Center,  and  the  corrected  data  are  directly  available  from  GETECH. 

All  three  data  sets  have  been  merged  to  produce  the  Bouguer  gravity  map  of 
Northern  Africa  (Figure  2).  Details  on  the  availability  of  these  data  can  be  found  in 
Appendix  B. 

Seismic  Data  Base 

We  have  selected  a  series  of  magnitude  >  5.0  earthquakes  and  their  associated 
aftershocks  occurring  between  1963  and  1988  (the  major  years  of  operation  of  the 
WWSSN)  for  study  (Figure  3,  Table  1).  Many  of  the  mainshocks  have  focal 
mechanisms  determined  from  teleseismic  waveform  modeling  or  first  motion  studies 
(Table  2).  Some  sequences  produced  numerous  aftershocks,  and  we  are  collecting  data 
for  all  events  above  magnitude  3.5.  The  data  we  are  collecting  consist  of:  1 )  location, 
origin  time,  arrival  time,  delay  time,  magnitude,  and  first  motion  information  from 
listings  of  the  International  Seismological  Centre  or  the  U.S.  Geological  Survey's 
Earthquake  Data  Report,  2)  published  focal  mechanisms,  focal  depth  and  other  source 
parameter  information  determined  from  teleseismic  or  local  network  information,  and  3) 
waveform  data  copied  from  microfiche  for  WWSSN  stations  located  at  regional  distances 
from  the  events  of  interest. 

We  have  collected  about  75%  of  the  listing  information  and  have  entered  much  of 
the  information  into  a  PC  based  version  of  the  EXCEL  spreadsheet/data  base  software. 
Examples  of  spreadsheets  and  other  details  related  to  this  information  are  given  in 
Appendix  C.  Figure  4  shows  an  example  of  how  we  plan  to  compare  P-wave  delay 
times  versus  azimuth  for  events  within  different  regions  of  North  Africa.  This  plot 
shows  the  delays  for  a  magnitude  4.5  aftershock  of  the  1969  northern  Red  Sea 
earthquake  sequence. 

We  have  conducted  a  fairly  extensive  search  of  the  literature  for  focal  mechanism 
information  and  believe  we  have  collected  about  90%  of  the  available  focal  mechanisms. 
We  hope  to  determine  focal  mechanisms  for  other  M>  5.0  events  from  first  motion 
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Bouguer  Gravity  Map  of  Northern  Africa 
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Figure  2.  Preliminary  Bouguer  gravity  map  of  North  Africa.  Contour  interval:  20mgal.  Precambrian  outline  of  Ahaggar  (Algeria), 
Tibesti  (S.  Lybia),  and  Darfur  (E.  Chad  and  W.  Sudan)  domal  uplifts  shown  in  dark  black.  Locations  of  three  crustal- scale  gravity  ’ 


Figure  3.  Magnitude  >  5.0  earthquakes  of  North  Africa  (1963-1988).  Large  circles  are 
magnitude  >  6.0  earthquakes.  Labeled  diamonds  represent  locations  of  WWSSN 
stations  operating  during  the  period  of  interest. 
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Table  1.  Magnitude  >  5.0  Earthquakes  (1963-1988)  in  North  Africa  of  Interest  to  This 
Study 


Date 

Time 

Magnitude 

Location 

Feb.  21,  1963 

1724 

5.0 

Gulf  of  Sidra  1 

May  11,  1963 

0110 

5.8 

Gulf  of  Sidra 

Jan.  1,  1965 

2138 

5.0 

M'Sila,  Algeria 

Jan.  26,  1967 

1611 

5.1 

Tunisia 

July  13,  1967 

0210: 

5.0 

Algeria 

Apr.  17,  1968 

0912 

5.1 

Boudinar,  Morocco 

Mar.  24,  1969 

1154 

5.2 

N.  Red  Sea,  Egypt 

Mar.  31,  1969 

0715 

6.0 

N.  Red  Sea,  Egypt 

Mar.  31,  1969 

2142 

5.0 

N.  Red  Sea,  Egypt 

Apr.  8,  1969 

1031 

5.2 

N.  Red  Sea,  Egypt 

Apr.  16,  1969 

0812 

5.0 

N.  Red  Sea,  Egypt 

Apr.  23,  1969 

1337 

5.0 

N.  Red  Sea,  Egypt 

Dec.  1,  1970 

0102 

5.2 

Tunisia 

Jan.  12,  1970 

0815 

5.1 

N.  Red  Sea,  Egypt 

Sept.  2,  1972 

1453 

5.3 

Gulf  of  Sidra 

June  28,  1972 

0949 

5.5 

N.  Red  Sea,  Egypt 

Nov.  24,  1973 

1405 

5.1 

Algeria 

Nov.  24,  1973 

1522 

5.3 

Algeria 

April  29,  1974 

2004 

5.0 

Nile  Delta,  Egypt 

Sept.  4,  1974 

0629 

5.5 

Tripoli,  Libya 

Jan.  1,  1977 

2046 

5.7 

A'in  Soltane,  Tunisia 

Dec.  9,  1978 

0713 

5.3 

W.  Egypt 

Oct.  10,  1980 

1539 

6.1 

El  Asnam,  Algeria 

Oct.  10,  1980 

1733 

5.2 

El  Asnam,  Algeria 

Oct.  10,  1980 

1908 

5.0 

El  Asnam,  Algeria 

Oct.  13,  1980 

0637 

5.0 

El  Asnam,  Algeria 

Oct.  14,  1980 

1734 

5.0 

El  Asnam,  Algeria 

Oct.  30,  1980 

2338 

5.1 

El  Asnam,  Algeria 

Nov.  8,  1980 

0754 

5.1 

El  Asnam,  Algeria 

Dec.  12,  1980 

1332 

4.8 

El  Asnam,  Algeria 

Dec.  7,  1980 

1737 

5.1 

El  Asnam,  Algeria 

Jan.  15,  1981 

0425 

5.1 

El  Asnam,  Algeria 

Feb.  1,  1981 

1320 

5.2 

El  Asnam,  Algeria 

Oct.  27,  1985 

1934 

5.7 

Constantine,  Algeria 
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Table  2.  Focal  Mechanisms  for  M>5.0  North  African  Earthquakes  (1963-1988) 


Date 

Mechanism 
strike. dip. rake 

Data 

Used 

Reference 

Feb.  21,  1963 

243,75,90 

first  motion 

McKenzie  (1972) 

Jan.  1,  1965 

70,65,68 

first  motion 

McKenzie  (1972) 

July  13,  1967 

167,30,90 

first  motion 

McKenzie  (1972) 

Mar.  31,  1969 

330,50,5 

waveforms 

Ben-Menahem  et  al.  (1976) 

325,50,-40 

first  motion 

McKenzie  et  al.  (1970) 

Jan.  12,  1972 

8,86,-10 

waveforms 

Pearce  (1980) 

June  28,  1972 

350,80,10 

waveforms 

Pearce  (1977) 

Apr.  29,  1974 

310,47,-17 

waveforms 

Ben-Menahem  et  al.  (1976) 

Sept.  4,  1974 

297,37,-141 

waveforms 

Westaway  (1990) 

Oct.  10,  1980 

225,54,83 

waveforms 

Deschamps  et  al.  (1982) 

Oct.  30,  1980 

42,40,72 

first  motion 

Yielding  et  al.  (1989) 

Nov.  8,  1980 

41,64,75 

first  motion 

Yielding  et  al.  (1989) 

50,30,97 

first  motion 

Ouyed  et  al.  (1983) 
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information  we  have  collected  using  the  grid  search  program  of  Whitcomb  (1973)  to 
determine  a  best  fit  focal  mechanism. 

About  50%  of  the  regional  waveform  data  have  been  collected,  and  we  have  just 
begun  to  digitize  these  seismograms.  Figure  5  shows  an  example  of  an  M=5.0 
aftershock  of  the  October  1980  El  Asnam  sequence  recorded  at  Toledo,  Spain  (TOL) 
located  at  625  from  the  event.  During  the  second  year  of  the  research  we  will  model 
these  waveforms  to  validate  and  update  the  preliminary  lithospheric  models  obtained 
from  the  gravity  data.  We  will  employ  reflectivity  (Randall  and  Taylor,  written  comm., 
1993)  and  full  waveform  finite  difference  (Sandmeier,  1990)  in  this  modeling  process. 


Results  and  Discussion 

Gravity  Modeling 

Once  we  had  merged  the  GETECH  data  into  our  existing  database,  we  selected 
three  profiles  across  the  region  to  model  (Figure  2).  These  profiles  represent  important 
propagation  paths  between  earthquakes  and  seismograph  stations  located  within  North 
Africa.  Profile  1  extends  from  central  Egypt,  the  site  of  several  recent  earthquakes  and 
seismograph  stations  (such  as  HLW),  across  central  Libya,  where  we  have  good 
geological  and  geophysical  control  on  crustal  structure  (e.g.  Suleiman,  1993;  Nyblade  et 
al.,  1996),  to  the  Hoggar  Uplift  of  southern  Algeria  (site  of  station  TAM,  Tamanrasset, 
and  several  French  nuclear  tests  in  the  early  1960's).  Profile  2  extends  across  the  Atlas 
Mountains  from  the  Atlantic  coast  of  Morocco  to  central  Tunisia.  Seismograph  stations 
are  located  along  this  profile  in  Morocco,  Algeria  and  Tunisia,  and  earthquakes  occur 
throughout  the  region.  Profile  3  extends  from  north  central  Algeria  across  the  Atlas 
Mountains  to  the  Hoggar  uplift.  This  profile  is  along  the  propagation  path  between 
northern  Algerian  earthquakes  and  TAM.  A  preliminary  velocity  model  has  been 
obtained  by  McNamara  et  al.  (1996)  from  regional  seismograms  whose  waves  propagate 
along  paths  that  parallel  Profile  3. 

The  preliminary  density  model  we  have  obtained  for  Profile  1  is  shown  in  Figure 
6.  The  upper  part  of  our  model  reflects  the  basin  structure  of  Egypt,  Libya  and  Algeria 
obtained  from  published  geological  and  geophysical  studies  of  these  basins  (e.g. 
Goudarzi,  1980;  van  der  Meer  and  Cloetingh,  1993;  Guirard  and  Maurin,  1992;  Exxon, 
1985),  primarily  related  to  their  petroleum  potential.  Since  our  knowledge  of  the  upper 
crust  below  the  sedimentary  basins  and  the  lower  crust  is  poor,  we  have  treated  these 
units  as  having  uniform  densities.  We  varied  thicknesses  of  these  units  only  as  required 
to  match  the  long  wavelength  features  of  the  observed  gravity  data.  Note  that  we  have 
not  attempted  to  match  every  short  wavelength  feature  of  the  profile,  since  our 
knowledge  of  the  shallow  geology  is  incomplete.  The  shallow  geology  also  should 
influence  regional  waveform  propagation  less  than  the  deeper  crustal  structure. 

Profile  1  shows  a  relatively  uniform  boundary  between  the  upper  and  lower 
crust.  There  is  a  gradual  shallowing  of  the  crust/mantle  boundary  between  the  Hoggar 
Uplift  (38  km)  and  the  Egyptian  coast  of  the  Mediterranean  (26  km)  (Sweeney,  1995). 

A  similar  eastward  shallowing  of  the  crust/mantle  boundary  is  observed  in  northern 
Libya  (Suleiman,  1993).  The  crustal  thickness  of  the  Hoggar  Uplift  is  in  good 
agreement  with  the  velocity  model  of  McNamara  et  al.  (1996)  and  recent  receiver 
function  studies  at  TAM  (E.  Sandoval,  pers.  comm.,  1996). 
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Figure  6.  Preliminary  interpretation  of  gravity  data  along  Profile  1  trending  West-East 
from  Ahaggar  domal  uplift  (Algeria)  to  Mediterranean  coastline  (Egypt)  (total 
length=2620  km).  Above:  Observed  and  calculated  gravity  values  for  each  station  along 
profile.  Below:  Crustal  model  used  to  generate  calculated  gravity  values.  Moho 
discontinuity  shown  in  bold.  Numbers  are  densities  for  the  bodies  (g/cm^). 
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Figure  7  shows  a  preliminary  interpretation  of  Profile  3.  The  shallow  upper 
crustal  structure  is  constrained  from  knowledge  of  basin  geometries  and  thicknesses 
(e.g.  Schurmann,  1974;  UNESCO,  1971;  Lesquer  et  al.,  1988).  The  Moho  depth  at  the 
northern  end  of  the  profile  is  constrained  by  the  results  of  the  European  Geotraverse 
(Ansorge  et  al.,  1992).  The  crust  appears  to  rapidly  thicken  beneath  the  Atlas 
Mountains,  then  gradually  thicken  to  the  south.  Again,  the  Hoggar  Uplift  can  be 
modeled  with  a  38  km  thick  crust.  A  gravity  low  over  the  uplift  suggests  the  presence  of. 
an  upper  mantle  anomaly  with  width  of  200  to  300  km  and  depth  of  40  to  60  km.  This  is 
much  smaller  than  observed  under  the  Ethiopian  and  Kenya  domes. 

Crustal  Thickness  Map 

Figure  8  shows  a  preliminary  map  of  crustal  thickness  based  on  an  integrated 
analysis  of  the  limited  seismic  constraints  and  Bouguer  gravity  data.  Note  that  the  crust 
is  generally  35-40  km  thick  away  from  the  continental  margins.  There  is  a  moderate 
amount  of  crustal  thickening  associated  with  the  Atlas  Mountains.  Crustal  thinning  is 
associated  with  the  Sirt  Basin  and  the  West  and  Central  African  rift  system  and  the 
adjecent  Cameroon  line.  The  lack  of  a  major  crustal  anomaly  associated  with  the  domal 
uplifts  such  as  the  Hoggar  is  suprising,  and  we  will  be  investigating  this  and  updating 
the  contour  map  as  our  study  progresses. 
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Figure  7.  Preliminary  interpretation  of  gravity  data  along  Profile  3  trending  South-North 
from  Air  massif  (Niger)  to  Mediterranean  coastline  (Algeria)  (total  length -224U  km). 
Above:  as  in  Figure  6.  Below:  as  in  Figure  6,  with  body  of  density  3.25  g/cm 
representing  possible  mantle  upwelling. 
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Conclusions 


We  are  continuing  our  modeling  of  gravity  Profile  2  along  the  axis  of  the  Atlas 
Mountains.  Further  2-D  filtering  and  analysis  of  the  gravity  data  set  is  also  planned. 
Modeling  of  regional  seismic  phases  will  continue  throughout  the  next  year.  We  will 
analyze  the  variations  in  travel  times  and  magnitude  estimates  for  the  earthquake 
sequences  of  interest.  Our  final  results  will  be  a  geophysical  and  geological  data  base  for 
North  Africa  that  will  be  easily  accessible  to  the  scientific  community  and  a  series  of 
"best”  lithospheric  models  for  regions  of  North  Africa  that  have  been  validated  through 
waveform  modeling  of  regional  seismic  phases. 
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Appendix  A  -  Documentation  of  Geological  Data  Base 

The  following  pages  show  maps  of  North  African  geological  features  that  we 
have  digitized.  The  data  are  available  at  the  UTEP  Department  of  Geological  Sciences 
web  site  under  North  Africa  Project.  The  sources  of  data  used  for  the  maps  we  have 
generated  are  given  below,  along  with  information  on  the  scales  of  the  original  maps. 

Figure  2  shows  the  outlines  of  exposed  Precambrian  rocks  in  the  Hoggar,  Tibesti 
and  Darfur  uplifts  (heavy,  dark  lines).  These  exposures  were  taken  from  UNESCO 
(1986)  (1:5,000,000). 

Two  sets  of  data  are  shown  in  Figure  Al.  The  white  lines  in  the  south  and 
northeastern  portions  of  the  study  area  represent  Early  Cretaceous  basins  and  faults 
taken  from  Guirard  and  Maurin  (1992)  (-1:55,000,000).  White  contours  in  northern 
Tunisia  show  a  Moho  contour  map  from  Research  Group  for  Lithospheric  Structure  in 
Tunisia  (1992)  (-1:600,000). 

Rift  basins  in  Niger,  Chad  and  the  central  African  Republic  are  shown  by  white 
lines  in  Figure  A2.  They  have  been  digitized  from  Genik  (1992)  (-1 :2,500,000). 
Basement  faults  in  north-central  Libya  are  from  van  der  Meer  and  Cloetingh  (1993). 

Figure  A3  shows  structural  features  of  central  Libya  from  Goudarzi  (1980) 
(-1:1,500,000).  Structure  in  Algeria  and  Chad  is  from  Dautria  and  Lesquer  (1989) 
(1:10,500,000). 

Locations  and  thickness  of  Early  Cretaceous  basins  in  the  northeastern  portion  of 
North  Africa  is  shown  on  Figure  A4.  These  data  are  from  Van  Houten  (1983) 
(-1:1,500,000). 

Contour  map  at  base  of  Mesozoic  rocks  in  Libya  (Goudarzi,  1980; 

-1:1,500,000)  is  shown  in  Figure  A5.  Major  geotectonic  units  of  West  Africa  (Lesquer 
et  al.,  1990)  are  also  shown. 

A  structural  contour  map  of  the  top  of  the  Precambrian  basement  in  Libya 
(Goudarzi,  1980;  -1,500,000)  is  shown  in  Figure  A6.  Lines  in  western  portion  of 
North  Africa  are  structures  shown  in  Kirton  and  Bourennani  (1984)  (-1:3,000,000). 
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ouguer  Gravity  Map  of  Northern  Africa 
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Figure  Al,  Cretaceous  basins  and  faults,  Moho  contours  for  Tunisia 


Bouguer  Gravity  Map  of  Northern  Africa 
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Figure  A2.  Rifts  of  southern  study  area,  structure  of  Libya 


Bouguer  Gravity  Map  of  Northern  Africa 
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Figure  A3.  Structure  of  central  Libya,  structure  of  Hoggar  uplift 


Bouguer  Gravity  Map  of  Northern  Africa 


Figure  A4.  Cretaceous  basins  of  northeastern  study  area 


ouguer  Gravity  Map  of  Northern  Africa 
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Figure  A5.  Structure  of  Mesozoic  in  Libya,  tectonic  features  of  West  Africa 


Bouguer  Gravity  Map  of  Northern  Africa 
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Figure  A6.  Structure  of  Precambrian  in  Libya,  structure  of  West  Africa 


Appendix  B  -  Documentation  of  Gravity  Data  Base 

Two  gravity  data  bases  have  been  compiled  for  public  use  in  North  Africa.  The 
first  of  these  in  the  Defense  Mapping  Agency  (DMA)  data  base  that  includes  36,410 
gravity  points  distributed  over  North  Africa  north  of  10  degrees  north.  These  data 
include  principal  facts  and  are  in  the  following  format:  field  #1  -  latitude  in  decimal 
degrees,  field  #2  -  longitude  in  decimal  degrees,  field  #3  -  Bouguer  anomaly  in  milligals 
using  a  2.67  g/cc  reduction  density,  field  #4  -  inner  zone  terrain  correction.  Hammer 
zones  (A-F),  in  milligals,  field  #5  -  outer  zone  terrain  correction,  Hammer  zones  G  and 
beyond,  in  milligals,  field  #6  -  elevation  above  sea  level  in  feet,  field  #7  -  observed 
gravity  in  milligals,  and  field  #8  -  source  number  indicating  the  origin  of  the  data  point. 

A  value  of  -99.99  milligals  in  the  terrain  correction  fields  means  that  the  terrain  correction 
for  that  point  has  not  been  made.  This  data  base  is  available  on  the  UTEP  web  page 
under  "DMA  gravity". 

The  second  data  base  available  is  more  detailed  data  for  an  area  surrounding  the 
Gulf  of  Sirt.  Except  for  92  points  in  this  data  set  principal  facts  are  not  available  for 
these  data.  This  data  set  includes  5531  gravity  points.  The  format  is  the  same  as  for  the 
DMA  gravity  data  described  above,  except  that  only  the  first  three  fields  of  data  are 
given.  This  data  base  is  available  on  the  UTEP  web  page  under  "Sirt  gravity". 

A  third  proprietary  data  base  was  used  to  prepare  the  maps  in  this  report  and 
those  available  on  the  UTEP  web  page.  This  data  base  contains  477,298  gravity  points, 
including  those  virtually  all  of  the  data  points  in  the  two  data  sets  mentioned  above  and 
has  been  gridded  at  5  minute  intervals.  Neither  the  data  points  nor  the  gridded  data  are 
available  to  the  public.  Individuals  with  specific  interest  in  these  data  should  contact  G. 
Randy  Keller  at  UTEP. 
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Appendix  C  -  Documentation  of  Seismic  Data  Base 

We  have  collected  earthquake  listings  from  the  bulletins  of  the  International 
Seismological  Centre  and  the  U.S.  Geological  Survey’s  Earthquake  Data  Reports  for 
earthquakes  of  M>5.0  and  associated  aftershocks  occurring  between  1963  and  1988 
within  North  Africa.  These  data  have  been  organized  into  a  data  base  using  the  PC  based 
program  EXCEL.  We  have  chosen  this  format  for  two  reasons.  First,  it  allows  our 
students  the  ability  to  work  at  home  at  odd  hours  to  enter  information  in  the  data  base. 
Second,  the  students  can  relatively  easily  make  relational  plots  from  the  data  base.  Two 
diskettes  containing  the  data  base  have  been  included  with  the  draft  of  this  report. 
Additional  diskettes  are  available  from  us  upon  request.  We  eventually  plan  to  load  the 
data  base  onto  our  departmental  main  frame  so  that  it  can  be  accessed  by  users  through 
the  Internet. 

Figure  Cl  shows  an  example  of  the  main  spread  sheet  we  have  created  for  all 
events  of  interest.  The  spread  sheet  contains  information  on  the  date,  time,  latitude  (lat), 
longitude  (Ion),  depth,  magnitude  (Ms,  mb  or  other),  and  location.  We  also  keep  track 
of  whether  we  have  obtained  listings  for  the  earthquake  and  plan  to  add  columns  for 
focal  mechanism  information. 

For  each  earthquake  for  which  we  have  obtained  a  listing,  we  have  created  the 
spreadsheet  shown  in  Figure  C2.  The  first  line  of  the  spreadsheet  contains  header 
information  on  the  earthquake  itself.  Then  we  list  all  stations  within  20°  of  the 
earthquake  and  note  the  station  distance,  station  time,  P-wave  arrival  time,  S-wave 
arrival  time,  and  arrival  time  of  other  phases,  if  given.  We  also  record  the  delay  time 
(observed  arrival  time  minus  predicted  time).  Magnitude  information  (usually  log  A/T)  is 
also  noted,  when  available.  We  will  add  another  column  that  will  give  information  on 
whether  we  were  able  to  obtain  seismograms  for  a  particular  station  and  the  status  of  the 
seismograms  (e.g.  digitized,  processed,  modeled)  when  they  have  been  obtained. 


27 


LISTING 

X 

X 

X 

X 

X 

X 

1 

X 

X 

i 

1  X 

i 

i 

X 

X 

1 

X 

X 

X 

X 

X 

X 

! 

X 

X 

X 

i 

i* 

i * 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

j 

i  X 

i 

! 

* 

I 

1 

! 

<D 

C 

ro 

o 

CO 

cr 

< 

i 

i 

i 

r~ 

r- 

1 

i 

1 

i 

1 

; 

j 

j 

| 

i 

i 

i 

o 

1- 

< 

o 

o 

_J 

Q. 

>. 

O) 

LU 

Q. 

>. 

o> 

LU 

CL 

>, 

03 

LU 

CD 

i_ 

ID 

CO 

CD 

%— 

*o 

CO 

CD 

2 

CO 

CP 

c 

o 

CO 

c 

< 

iE 

CD 

C 

m 

< 

LU 

El  Asnam 

I 

c 

CO 

LS 

£ 
CD 
C 
c n 
< 

LU 

l 

IE 

CD 

c 

If) 

< 

,LU 

E 

CD 

C 

CO 

< 

LU 

E 

CD 

C 

CO 

< 

LU 

E 

< 

in 

E 

2 

CO 

< 

LU 

E 

CD 

c 

CO 

< 

LU 

E 

CD 

C 

CO 

< 

LU 

E 

CD 

C 

CO 

< 

LU 

1  E 

CD 

C 

CO 

< 

\UJ 

E 

CD 

c 

CO 

< 

LU 

£ 

CD 

C 

CO 

< 

UJ 

E 

2 

« 

5 

LU 

E 

CD 

c 

CO 

< 

ED 

E 

CD 

C 

CO 

< 

LU 

E 

CD 

C 

CO 

< 

Hi 

E 

CD 

C 

CO 

< 

LU 

E 

CD 

C 

CO 

< 

E 

CD 

C= 

CO 

< 

LU 

E 

CD 

C 

CO 

< 

S3 

E 

ca 

5 

< 

UJ 

E 

CD 

C 

CO 

< 

LU 

CD 

'i _ 

<D 

O' 

< 

If 

!S 
!< 
j  UJ 

111 

Q 

3 

H* 

2 

O 

< 

s 

-Q 

<fi 

* 

O 

in 

CD 

xT 

o 

lo 

CO 

in 

in 

in 

LO 

xr 

LO 

CD 

xT 

ID 

lO 

CD 

xr 

CO 

to 

5 

<q 

XT 

03 

X — 

CO 

o 

xr 

CM 

LQ 

O 

LO 

CO 

xT 

CM 

CM 

CO 

xr 

CO 

CM 

* 

CD 

xr 

CO 

xr 

r- 

xT 

CM 

xr 

CM 

xr 

LD 

<9 

O 

xr 

xr 

q 

xr 

jo 

id 

1 

L 

5 

to 

LO 

XT 

CD 

xT 

1? 

CO 

CD 

CD 

LO 

xT 

^r 

CD 

00 

CO 

xr 

oo 

X 

CO 

CD 

CO 

LO 

xr 

LO 

xr 

* 

CM 

xr 

CD 

CO 

CD 

xr 

q 

xr 

CD 

CM 

* 

CM 

LO 

CD 

CO 

xr 

juo 

| 

CM 

r-~’ 

CD 

‘p 

i 

xr 

CO 

X 

" 

CD 

CO 

X 

xr 

« 

-J 

CM 

LO 

xr 

CO 

jcq 

1  xr 

i 

i 

i 

DEPTH 

o 

CO 

CO 

9 

CM 

xr 

18.0 

o 

CO 

o 

o 

10.0 

21.0 

q1 

xr 

o 

o 

o 

d 

o 

£ 

o 

o 

o 

o 

o 

o 

o 

° 

O 

CO 

o 

o 

o 

d 

o 

o 

o 

o 

O  iO 
0:0 

i 

o 

o 

O 

d 

o 

o 

^o1 

d 

x — 

o 

o 

o 

d 

O 

CM 

o 

d 

10.0 

10.0 

o 

o 

§ 

CD 

CO 

CO 

CO 

CO 

CO 

CO 

o 

CD 

q 

CO 

o> 

CO 

XT 

CO 

xr 

CD 

xr 

LO 

cq 

X — 

in 

CO 

£ 

CD 

■’c” 

CD 

CD 

in 

x — 

N*  CO 

X 

. 

in 

T — 

q 

LO 

x — 

q 

q 

T — 

q 

q 

xr 

x — 

xr 

CD 

X — 

£ 

_i 

I 

27.3 

uq 

CM 

LO 

CM 

31.3 

33.0 

o 

CO 

CO 

36.5 

36.4 

CD 

CO 

36.3 

CD 

CO 

36.3 

o 

CD 

CO 

o 

CD 

CO 

^.i 

CD 

CO 

i 

J 

36.2 

CD 

CO 

CD 

CO 

36.4 

36.5 

36.2 

36.2 

CM 

CD 

CO 

CO 

CD 

CO 

CO 

CD 

CO 

xT 

CD 

CO 

36.3 

CO 

CD 

CO 

CO 

CD 

CO 

36.3 

36.0 

36.3 

36.3 

TIME 

8:12:56 

8:01:04 

13:37:17 

14:53:38 

6:29:14 

21:21:52 

20:46:53 

13:55:00 

12:25:17 

12:37:08 

14:12:27 

14:44:51 

15:18:46 

15:25:36 

15:36:52 

15:39:06 

17:15:17 

O) 

LO 

d 

q 

X — 

19:07:59 

20:22:41 

21:21:18 

23:55:33 

23:59:58 

1:29:16 

5:40:48 

19:30 

21:26:13 

0:36:49 

4:02:15 

6:37:36 

14:33:37 

20:26:41 

1:07:39 

17:34:55 

DATE 

CD 

CD 

CD 

4/17/69 

4/23/69 

ZLIZI6 

9/4/74 

12/19/76 

1/19/77 

2/9/77 

10/10/80 

10/10/80 

10/10/80 

10/10/80 

O  i 
CO 

o 

51 

j 

10/10/80 

10/10/80 

10/10/80 

10/10/80 

o 

cq 

o 

o 

10/10/80 

10/10/80 

10/10/80 

10/10/80 

10/10/80 

10/11/80 

10/11/80 

I  10/11/80 

10/11/80 

10/12/80 

10/12/80 

10/13/80 

o 

cq 

CO 

o 

x — 

10/13/80 

o 

cq 

o 

10/14/80 

10/14/80 

28 


Figure  Cl.  Example  of  page  from  master  event  spreadsheet  that  contains  information  on  all  earthquakes  we  plan  to  study. 
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Figure  C2:  Example  of  spreadsheet  data  for  an  aftershock  of  the  1969  northern  Red  Sea  earthquake  sequence, 
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P.O.B.  2286,  NOLON  58122  ISRAEL 

MATTHEW  SIBOL 
ENSCO,  INC. 

445  PINEDA  COURT 
MELBOURNE,  FL  32940 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  1663,  MS  D460 
LOS  ALAMOS,  NM  87545 


LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

POBOX  808,  MS  L-195 
LIVERMORE,  CA  94551 
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JEFFRY  STEVENS 
MAXWELL  TECHNOLOGIES 
P.O.  BOX  23558 
SAN  DIEGO,  CA  92123 


BRIAN  SULLIVAN 
BOSTON  COLLEGE 
INSITUTE  FOR  SPACE  RESEARCH 
140  COMMONWEALTH  AVENUE 
CHESTNUT  HILL,  MA  02167 

NAFI TOKSOZ 

EARTH  RESOURCES  LABORATORY,  M.I.T. 
42  CARLTON  STREET,  E34-440 
CAMBRIDGE,  MA  02142 


FRANK  VERNON 

UNIVERSITY  OF  CALIFORNIA,  SAN  DIEGO 
SCRIPPS  INSTITUTION  OF  OCEANOGRAPHY  IGPP,  0225 
9500  GILMAN  DRIVE 
LA  JOLLA,  CA  92093-0225 

TERRY  WALLACE 
UNIVERSITY  OF  ARIZONA 
DEPARTMENT  OF  GEOSCIENCES 
BUILDING  #77 
TUCSON,  A Z  85721 

DANIEL  WEILL 
NSF 

EAR-785 

420 1  WILSON  BLVD.,  ROOM  785 
ARLINGTON,  VA  22230 

RU  SHAN  WU 

UNIVERSITY  OF  CALIFORNIA  SANTA  CRUZ 
EARTH  SCIENCES  DEPT. 

1156  HIGH  STREET 
SANTA  CRUZ,  CA  95064 

PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

POBOX  999,  MS  K5-12 
RICHLAND,  WA  99352 


JAMES  E.  ZOLLWEG 
BOISE  STATE  UNIVERSITY 
GEOSCIENCES  DEPT. 

1910  UNIVERSITY  DRIVE 
BOISE,  ID  83725 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
8725  JOHN  J.  KINGMAN  ROAD 
FT  BEL  VOIR,  VA  22060-6218  (2  COPIES) 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  1663,  MS  C335 
LOS  ALAMOS,  NM  87545 


DAVID  THOMAS 
I  SEE 

29100  AURORA  ROAD 
CLEVELAND,  OH  44139 


LAWRENCE  TURNBULL 

ACIS 

DCI/ACIS 

WASHINGTON,  DC  20505 


GREG  VAN  DER  VINK 
IRIS 

1616  N.  FORT  MEYER  DRIVE 
SUITE  1050 

ARLINGTON,  VA  22209 

LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  808,  MS  L-205 
LIVERMORE,  CA  94551 


JAMES  WHITCOMB 
NSF 

NSF/ISC  OPERATIONS/EAR-785 
4201  WILSON  BLVD.,  ROOM785 
ARLINGTON,  VA  22230 

JIAKANG  XIE 
COLUMBIA  UNIVERSITY 
LAMONT  DOHERTY  EARTH  OBSERVATORY 
ROUTE  9W 

PALISADES,  NY  10964 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  61 16 

MS  0750,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0750 

OFFICE  OF  THE  SECRETARY  OF  DEFENSE 
DDR&.E 

WASHINGTON,  DC  20330 


TACTEC 

BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 

COLUMBUS,  OH  43201  (FINAL  REPORT) 
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PHILLIPS  LABORATORY 
ATTN:  XPG 
29  RANDOLPH  ROAD 
HANSCOM  AFB,  MA  0 1 73 1  -3  0 1 0 


PHILLIPS  LABORATORY 
ATTN:  TSML 
5  WRIGHT  STREET 
HANSCOM  AFB,  MA  01731-3004 


PHILLIPS  LABORATORY 
ATTN:  GPE 
29  RANDOLPH  ROAD 
HANSCOM  AFB,  MA  01731-3010 


PHILLIPS  LABORATORY 
ATTN:  PL/SUL 
3550  ABERDEEN  AVE  SE 
KIRTLAND,  NM  871 17-5776  (2  COPIES) 
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